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Phase field models frequently provide insight to phase transitions, and are robust numerical 
tools to solve free boundary problems corresponding to the motion of interfaces. A body of prior 
literature suggests that interface motion via surface diffusion is the long-time, sharp interface limit 
of microscopic phase field models such as the Cahn-Hilliard equation with a degenerate mobility 
function. Contrary to this conventional wisdom, we show that the long-time behaviour of degenerate 
Cahn-Hilliard equation with a polynomial free energy undergoes coarsening, reflecting the presence 
of bulk diffusion, rather than pure surface diffusion. This reveals an important limitation of phase 
field models that are frequently used to model surface diffusion. 


A key problem in modelling phase transitions in ma¬ 
terials lies in linking macroscopic interfacial motion and 
mesoscopic dynamics. A common approach are phase 
field models, which replace a sharp interface with or¬ 
der parameters that are continuous across the interface. 
Phase field models can be constructed from systematic 
coarse graining of the microscopic Hamiltonian, and are 
often written as a gradient minimisation of certain micro¬ 
scopic free energy functional n As such, they provide 
the crucial link between microscopic interactions and the 
kinetics of phase separation and pattern-formation. 

Mesocopic phase field descriptions have been widely 
used in the literature as their numerical approximation is 
less complicated than the approximation of macroscopic 
descriptions based on sharp interfaces. By using a con¬ 
tinuous order parameter field, phase field approaches are 
versatile enough to capture topological changes, and re¬ 
place the numerically challenging task of interface track¬ 
ing with integration of a time-dependent partial differen¬ 
tial equation. Therefore, the phase field formalism is in¬ 
creasingly used as a numerical approximation for a wider 
class of free boundary problems than what the free en¬ 
ergy describes microscopically 

A particularly noteworthy class of free boundary prob¬ 
lem is when the velocity of the interface Vn is proportional 
to the surface Laplacian of the mean curvature k and 

Vn = MAsK, ( 1 ) 

where Ai is the mobility. Equation © is known as the 
surface diffusion flow, and has been used as a model for 
maw complex processes such as electromigration in met¬ 
als 0 , heteroepitaxial growth @ and more recently solid- 
solid dewetting 

The Cahn-Hilliard equation with degenerate mobility 
is the commonly used phase field model to approximate 
surface diffusion (e.g. |7H14|L where the order parameter 
u is conserved and satisfies (in dimensionless units) 

Mt = -V-j, j =-eM(u)VAi, (2a) 

e^i =—e^V'^u + f'{u), (2b) 

/(w) = ^(1 - M('u) = l-u/ (2c) 


where M(u) is the mobility function; j is the flux; /r is 
the chemical potential; e is the interfacial tension which 
determines the width of the interface, and f{u) is the bulk 
free energy. Throughout the paper, we will assume the 
no-flux condition n • j = 0, and the variational Neumann 
condition n ■ Vu = 0 on the boundaries of the solution 
domain. The two pure phases are denoted by u = ±1, 
and the interface is located at the contour m = 0. 

The precise form of the mobility is usually chosen on 
thermodynamic grounds [mil- For the lattice-gas en¬ 
tropy s{u) = ulogu -I- (1 — M)log(l — u), the Einstein 
relation stipulates that the mobility is related to the en¬ 
tropy function s{u) via M{u) = {d'^s/ = 1 — 
This motivates the choice of mobility in ([2]). 

We are interested in the long-time behaviour when the 
initial mixture has separated into regions where u is ei¬ 
ther close to 1 or to —1, except for regions of width 0(e) 
close to the interface over which u transitions between 
these two regions. Heuristically, the width of this inter¬ 
face layer decreases as e 0. One may think that if the 
mobility function is degenerate and vanishes at the pure 
phases, the flux normal to the interface is suppressed and 
therefore only surface diffusion via mass flux tangential 
to the interface can occur. However, this heuristic ar¬ 
gument neglects the fact that the gradient of the inter¬ 
face diverges as e —>■ 0. Therefore, whether the degener¬ 
ate mobility function is sufficient to suppress the normal 
mass flux at leading order is unclear [13, Hi ■ 

The key result of this paper is the presence of a nonlin¬ 
ear bulk diffusion term at leading order for the interface 
velocity, and the correct sharp interface limit that de¬ 
scribes the quasistationary evolution of the interface /, 
located at u = 0, is given by 


0, in H_, 

(3a) 

-At, on C, 

(3b) 

2 1 

-AsK + -^j.VnU on C, 

(3c) 

0, on 

(3d) 


where is the normal velocity, and the definitions of 
Ul- and C are given in Figure [TJ The possibility of bulk 
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FIG. 1: Illustration of the asymptotic structure of the degen¬ 
erate mobility case. The domain fiext is split into interfaces ^ 
where u = 0 which enclose the solid domains Q+ where u > 0 
(coloured red), and the region [2- with vapor (colored green) 
where u > 0. The normal n to the interfaces points out of 
II -I-. 


diffusion was noted in earlier works [i3[il and analysed 
recently by Dai and Dn [l^. However, in m the analysis 
was done on the (unphysical) solution branch with |u| > 1 
in some region. Our analysis below considers the physical 
branch of solution where |it| < 1 everywhere, and derives 
the limiting model as e —>■ 0. It can be shown using 
rigorous mathematical analysis that the physical branch 
with |n| < 1 everywhere exists for all parameter values 

M- 

To begin our analysis, we drop the time derivative from 
Rewriting the Laplacian in polar coordinates, and 
resolving the boundary layer near the interface by noting 
that the interface is 0(e) thick, we obtain the leading 
order solution 


( T — 7*0 \ I 

Mo(r)=tanh( —-— 1, /r =(4) 

Here r is the radial coordinate with respect to the centre 
of the osculating circle to the interface, tq is the posi¬ 
tion of the interface in these coordinates, and /i|^ is a 
constant. Equation (|3|) reveals how the interface is being 
represented as a continuous, albeit thin, order parameter 
profile of width e around r = tq. A key physical intuition 
is that the relaxation of the local order parameter profile 
to (31) is rapid (as it is driven by local rearrangement of 
particles), and the late stage dynamics is quasi-static and 
determined by the movement of the interface, i.e. change 
in tq. 

To put this intuition onto firmer footing, we consider 
the chemical potential ^ away from the interface. There, 
the system is almost a pure phase (say u = —1 + eu, 
/i = efl, with u, fi = 0(1)), and the spatial variation of 
u is negligible so that the Laplacian can be neglected. 
Expanding Equation (I2bl) in powers e, the leading order 


term is given by 

= (5) 

and substituting (O into (l2a)) . we obtain 

zii = M'(-1)V • (uV/i) = • (AV/i). (6) 

Further, we use the aforementioned quasi-static approx¬ 
imation, and assume that ut is small. Thus, to leading 
order. 


V • (AV/i) = 0. (7) 


The chemical potential at the interface can be obtained 
by multiplying (I2bl) by drU, and integrating both sides. 
Assuming that the azimuthal variations in u is asymp¬ 
totically smaller than the radial variation, we have 


rro+v r 

/ efidrudr = / 
J rn — v Jri 


ro+'T] - 


- dr{rdru) + f'{u) 


drU dr. 


( 8 ) 

Close to the interface, u can be approximated by u{r) 
uo{r). Substituting (3]) into dH]), and assuming e <C ?7 
1, we obtain 


I 2 

(9) 

where k = 1 /rg is the curvatnre of the interface. 

In the quasi-static approximation, we neglect time de¬ 
pendence except for the slow motion of the interface. To 
obtain the interface velocity, we focus on the boundary 
layer region, and move into a Lagrangian frame by mak¬ 
ing the transformation ut ^ Ut — {'v ■ V)u. Now, noting 
that normal velocity is much larger than the lateral ve¬ 
locity, we have (v • V)u « e^VndrU, where is the 
normal velocity, scaled in anticipation of the order of the 
right hand side which will turn out to be O(e^). Equation 
(El) becomes 

- e^VndrU=-dr{j-n) - ds{i-t), (10) 


where s denotes the coordinate tangent to the interface, 
and n and t are the (outward pointing) unit normal and 
the unit tangent vector, respectively. Close to the inter¬ 
face, the tangential flux is given by j • t = —eM{u)dsfi = 
—eM('u)9s^|^. Thus, again substituting (31) for u, and 
integrating over the interface, we arrive at 

2e^Vn = - [j ■ -b ^e^dssK. (11) 

Now, for e ^ ry <C 1, we obtain 


J ■ n 


ro+ri / (- 1 ) 


0(c3) 


and j • n|^ _^ = 0 as the system reaches the pure phase 
u = 1 deep in f2+. Therefore, all in all, we obtain (I5cl) . 
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e 

0.01 

0.005 

0.002 

0.001 

Eq (dlD 

Eq (HID 

a 

-133.2 

-133.8 

-136.3 

-137.0 

-137.4 

-128 


TABLE I: Relaxation rates for m = 2 obtained from the lin¬ 
earised phase field model are shown for different values of e in 
the first four columns, and compared to the eigenvalues ob¬ 
tained for linearised sharp interface models for pure surface 
diffusion m and the porous medium type model (1131) in the 
next-to-last and the last column, respectively. 


where we have identified dgs with the surface Laplacian 
Ag. Equation (15^ shows that the interface velocity in 
the sharp interface model has two contributions: one 
from surface diffusion, Equation o, which is local to 
the interface, and another contribution from nonlinear 
bulk diffusion which satisfies a porous-medium equation 
(I3a|) . Unlike pure surface diffusion, the mass flux aris¬ 
ing from bulk diffusion couples disjoint solid domains 
with each other. This results in coarsening where larger 
solid domains grow at the expense of smaller ones, which 
cannot happen for pure surface diffusion. Moreover, for 
non-circular interfaces, the contribution from surface and 
bulk diffusion enter the interface evolution to the same 
order, i.e. the effect of the latter does not become negli¬ 
gible compared to the former even when letting e —> 0. 

To test this prediction of our analysis, we consider 
the relaxation of an azimuthal perturbation to a radi¬ 
ally symmetric stationary state with radius tq and hence 
curvature k = I/tq. Eor azimuthal perturbations pro¬ 
portional to cos m6, the pure solid diffusion model © 
predicts an exponential decay rate 

a = —— 1 ). ( 12 ) 

3 

In contrast, the decay rate in the porous medium model 
© is given by 


(J = —K 


2 1 

—m?(inn? — 1) H—m(m^ — 1) tanh(m log k) 
3 9 


(13) 


Table U shows the decay rate numerically obtained by 
solving the phase field model. Equation ([5]). It shows 
how the decay rate of the azimuthal perturbation to the 
axisymmetric base state tends to the linearised sharp in¬ 
terface model with with the contribution from nonlinear 
bulk diffusion, rather than to the one for pure surface 
diffusion. 

Our analysis only applies to free energy functions / 
for which /'(±1) = 0 and /"(±1) > 0 at the minima 
M = ±1. Indeed, according to the asymptotic analysis in 
[^ . the double obstacle free energy 


/do(u) 


\ — V? if |m| < I, 
oo if |m| > I 


gives rise to pure surface diffusion flow for M{u) = I — 
V? in the sharp interface limit (e 0), and also if the 
logarithmic free energy f{u) = —v? + e°‘{u\ogu + (1 — 
u) log(I — u)) with a > 0 is used instead. 

In conclusion, our analysis establishes that phase field 
models for pure surface surface diffusion cannot be re¬ 
alised using the Cahn-Hilliard equation with the degen¬ 
erate mobility and Ginzburg-Landau free ener gy a s in ([3 , 
as was repeatedly assumed in the literature |9l.ll0l. ll3lll4l[ . 
A nonlinear bulk diffusion term appears to leading order 
of the sharp interface limit, hence affecting the coarsen¬ 
ing behaviour on the same time scale as surface diffusion. 
In particular, it allows, on this time scale, disjoint inter¬ 
faces to coarsen and cannot be suppressed by reducing 
e. 

We note that the derivation presented in this work 
could be made mathematically robust via matched 
asymptotic analysis. Such an approach was applied to 
analyse the Cahn-Hilliard equation with constant mobil¬ 
ity UM ■ Extending the method of matched asymp¬ 
totics to model © with a degenerate mobility will be the 
subject of a subsequent publication. The heuristic ap¬ 
proach presented here, however, reveals clearly the salient 
physics involved in the sharp interface limit. 
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